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Abstract 

We have studied the adsorption of gas molecules (CO, NO, NO2, O2, N2, CO2, and NH3) on 
graphene nanoribbons (GNRs) using first principles methods. The adsorption geometries, ad- 
sorption energies, charge transfer, and electronic band structures are obtained. We find that the 
electronic and transport properties of the GNR with armchair-shaped edges are sensitive to the 
adsorption of NH3 and the system exhibits n-type semiconducting behavior after NH3 adsorption. 
Other gas molecules have little effect on modifying the conductance of GNRs. Quantum transport 
calculations further indicate that NH3 molecules can be detected out of these gas molecules by 
GNR based sensor. 
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I. INTRODUCTION 



Sensing gas molecules is critical to environmental monitoring, control of chemical pro- 
cesses, space missions, and agricultural and medical applications 11. Solid-state gas sensors 

nn 

are renowned for their high sensitivity which have made them ubiquitous in the world[2|, |3|. 
In the past few years, a new generation of gas sensors have been demonstrated using carbon 
nanotubes (CNTs) and semiconductor nanowiresj^, [s, y, Q, 1^, lo, Q, 11]. It was reported 
that semiconducting CNTs could be used to detect small concentration of NH3, NO2, and 
O2 with high sensitivity by measuring changes of the CNTs conductance upon exposure to 



the gases at room temperature 



m 



GraphenCj a single atomic layer of graphite, has been successfully produced in 



experiments [Id , llSj , which have resulted in intensive investigations on graphene-based struc- 
tures because of fundamental physics interests and promising applications [l^ (e.g'., gas 
sensor 1^). Importantly, graphene can be patterned via standard lithographic techniques 



into quasi-one-dimension materials 13l. lid . llTj . graphene nanoribbons (GNRs), which have 



many properties similar to CNTs, such as energy gap dependence of widths and crystallo- 



graphic orientation [isl. Il9l| . Different from CNTs, however, GNRs present 



edges which make GNRs not only notably more accessible to doping 



18 



M 



on g and reactive 
and chemical 



modification 



mm 



2l|. 



;5| , but also more susceptible to structural defects and impurites 
Theoretical studies of gas molecular adsorption on the graphene surface have been re- 
ported recently 26|, l27|], which showed that NO2, H2O, NH3, CO and NO molecules are 
physically adsorbed on the pristine graphene: NH3 and CO molecules will act as donors 
while H2O and NO2 will act as acceptors, which are consistent with previous experiment jisl . 
Compared with graphene, GNRs are advantageous in small volume and free reactive edges. 
In experiments, the edges of GNRs are not well controlled 16|, ll7|] and it is hard to obtain 



fully saturated edges without any dangling bond (DB) defects. It is well known that DB 
defects around the vacancy sites or at the tips pla y a very important role in CNTs gas 
sensors because they are very chemically reactive |28l. Similar to CNTs, when there 

are DB defects at GNRs edges, covalent attachment of chemical groups and molecules also 
significantly influences their electronic properties 



2J, |25|. Therefore, it is very interesting 
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and important to study the feasibility of using GNRs as gas sensors. 

In this article, using density functional theory (DFT) calculations, we study the adsorp- 
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tion of gas molecules (CO, NO, NO2, O2, N2, CO2 and NH3) around the sites of DB defects 
on armchair GNRs (AGNRs, having armchair-shaped edges) and explore the feasibility of 



using AGNRs as gas sensors. Following conventional notation[31[, a GNR is specified by the 
numbers {n) of dimer lines and zigzag chains along the ribbon forming the width, for the 
AGNR and zigzag GNRs (ZGNRs, having zigzag-shaped edges), respectively. For example, 
the structure in Fig. la is referred as a 10- AGNR (i.e., n = 10). Previous works show that 



all AGNRs are semiconductor while ZGNRs are metal|18l. \3l\. We focus on semiconducting 
AGNRs instead of metallic ZGNRs, since it is expected that gas molecule adsorption will 
have a much smaller effect on modifying the electronic properties of (metallic) ZGNRs. It 
is found that although all gas molecules can influence the electronic structure of AGNRs, 
only NH3 molecule adsorption can modify the conductance of AGNRs remarkably by acting 
as donors while other gas molecules have little effect on conductance. This property can 
be utilized to detect NH3 out of other common gases, which is requisite and significant in 



industrial, medical, and living environments 321 ]. 



II. CALCULATION METHOD AND MODEL 

Our electronic structure calculations were performed using the DFT in the spin-polarized 
generalized gradient approximation (GGA) with PW91 functional for the exchange and 
correlation effects of the electrons, as implemented in Vienna Ab initio Simulation Package 



( VASP) [33| . The electron- ion interaction was described by the ultrasoft pseudopentials and 
the energy cutoff was set to be 400 eV. Structural optimization was carried out on all systems 
until the residual forces on all ions were converged to 0.01 eV/A. The quantum transport 



calculations were performed using the ATK 2.3 package 34| . which implements DFT-based 
real-space, nonequilibrium Green's function formalism. The mesh cutoff is chosen as 150 Ry 
to achieve the balance between calculation efficiency and accuracy. 

III. RESULTS AND DISCUSSION 

DB defects would exist at both edges of GNRs due to the fact that the two edges of a 
GNR are equivalent in nature. So we will focus our study on double-edge-defect case (i e., 
both edges have DB defects). What is more, we also preform test calculations for the single- 
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FIG. 1: (a) The optimized structure of 10-AGNR with edge danghng bond (DB) defects, where 
the arrow shows the periodic direction. The sites of DB defects are shown by yehow. C atoms and 
H atoms are denoted by large and small (white) spheres, respectively, (b) Density of states (DOS) 
of perfect 10-AGNR (top panel) and spin-polarized DOS of 10-AGNR with DB defects (middle 
and bottom panels). The Fermi level is set to zero (the top of valence band). The middle (blue) 
and bottom (red) panels in the figure correspond to minority spin and majority spin, respectively. 
The green area corresponds to local DOS of the two carbon atoms (yellow balls in (a)) with DB 
defects. 



edge-defect case (DB defects only exist at one edge), and find little difference in essential 
results compared with the double-edge-defect case. The structure with one DB defect per 
edge per five unit cells in the ribbon axis direction is adopted in our calculations, as shown 
in Fig. la, corresponding to a DB defect concentration of 0.04 A~^. [Several different initial 
configurations with one DB defect per edge have been considered and we find that the one 
shown in Fig. la is most favorable (stable) in energy.] Due to the dangling a-bonds at the 
edges, the ground state of this system is spin-polarized with the magnetic moment of 1 /xB 
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FIG. 2: Optimized structure of 10-AGNRs with gas molecule adsorption: (a) CO, (b) NO, (c) 
NO2, (d) O2, (e) CO2, and (f) NH3. We only show the structure around the adsorbed molecule. 

per DB, which is localized at the carbon atoms with DB defects. Spin-polarized density of 
states (DOS) of this system is shown in the middle and bottom panels of Fig. lb. For better 
comparison, the DOS of prefect 10-AGNR with the same supercell is also presented (the 
top panel). It can be seen that the perfect 10-AGNR is paramagnetic semiconductor with 
a band gap of 1.2 eV, consistent with our previous study ISj. When there are DB defects 
at the ribbon edges, two new peaks appear in the DOS (indicated by the arrows in Fig. 
lb). The local density of states (LDOS) analysis (the green area in Fig. lb) shows that the 
two peaks are mainly contributed by the edge carbon atoms with DB defects. These two 
DB states are localized, respectively, within the valence band for minority spin and at the 
bottom of conduction band for majority spin, different from the usual DB states which are 
localized within band gap 3^. 

We start by investigating the adsorption geometries of seven gas molecules on 10-AGNR 
with DB defects (Fig. la). Only the gas molecules adsorption around the sites of DB 
defects is considered. Several different initial orientations of gas molecules on 10-AGNRs 
are adopted in searching the most stable configurations. Fig. 2 shows the top views of 
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TABLE I: Calculated adsorption energies (-Eads) and charge transfer (CT) from the gas molecules 



to the 10-AGNR. 



Molecules 


CO 


NO 


NO2 


O2 


N2 


CO2 


NH3 


E^ds (eV) 


-1.34 


-2.29 


-2.70 


-1.88 


0.24 


-0.31 


-0.18 


CT{e) 


-0.30 


-0.55 


-0.53 


-0.78 


/ 


-0.41 


0.27 



10-AGNR with adsorbed molecules. We find that different gas molecules prefer different 
geometries in the adsorption: (a) CO molecule lies in the AGNR plane with C-C distance 
of 1.35 A and a C-C-0 angle of 168.6°. The bond length of the adsorbed CO is 1.18 A, a 
little longer than that of an isolated molecule (1.13 A), indicating that adsorption process 
will weaken the original C-0 bond of gas molecule, (b) NO molecule sits out of the AGNR 
plane with a C-N-0 angle of 118.2°, and the C-N and N-0 distances are 1.43 A and 1.24 
A, respectively. Also, the bond length of the adsorbed NO is 0.07 A longer than that of 
an isolated molecule, (c) The adsorbed NO2 has C-N and N-0 distances of 1.47 A and 
1.25 A and a 0-N-O angle of 127°, and the 0-N-O plane is tilted 60° with respect to the 
ribbon plane. Our results are consistent with previous studies 2^, l25|. (d) O2 sits in the 
ribbon plane with C-0 and 0-0 distances of 1.37 A and 1.38 A respectively and a C-0-0 
angle of 116.2°. The bond length of O2 molecule increases by 0.15 A after adsorption, 
(e) The geometry of CO2 adsorbed on the ribbon edge is different from that of an isolated 
CO2 molecule: it is not linear structure anymore but similar to the NO2 configuration. This 
indicates that the bond of CO2 transfers from sp into sp^ hybridization in order to lower 
the total energy. The C-C and C-0 distances are 1.51 A and 1.26 A respectively, and 
the C-O-C angle is 127°. (f) The adsorbed NH3 molecule sits 1.49 A away from the edge 
carbon atom, the H-N distance is 1.03 A and the dihedral angle is about 21° between the 
C-N bond and the AGNR plane. Besides the above gas molecules, we also study the N2 
molecule adsorption on the AGNR, and find that it is difficult for N2 to adsorb on the ribbon 
due to its inert nature. 

The calculated adsorption energies of the gas molecules with the 10-AGNRs are shown in 
Table 1. Herein, The adsorption energy is defined as: i^ads = {-^tot (ribbon + mMolecule) — 
i?tot (ribbon) — nE'tot (Molecule) }/m, where E'tot (ribbon + mMolecule), -Etot (ribbon), and 
i^tot (Molecule) are the total energies of the AGNR with molecule adsorption, the isolated 
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FIG. 3: Band structure and density of states (DOS) of 10-AGNRs with gas molecule adsorption: 
(a) CO, (b) NO, (c) NO2, (d) O2, (e) CO2, and (f) NH3. The LDOS of gas molecules is also plotted 
(red filled area under DOS curve). The Fermi level is set to zero. 

AGNR (with DB defects) and the molecules, respectively. And m is the number of molecules 
adsorption on the AGNR. The results reveal that all these adsorption configurations are ener- 
getically favorable except that the N2 adsorption process is endothermic reaction (note: the 
negative adsorption energy corresponds to the exothermic reaction). The adsorption energies 
of CO, NO, NO2, and O2 are all larger than 1 eV, corresponding to strong chemisorption. The 
adsorption energies for the NH3 and CO2 on AGNRs are -0.18 eV and -0.31 eV respectively, 
indicating that the adsorption are between weak chemisorption and strong physisorption. 
The above results illuminate that gas molecule adsorption at AGNR edges is quite different 
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FIG. 4: Iso-surface plots of the partial charge density at T point of the band crossing the Fermi 
level for the 10-AGNRs with (a) CO2 adsorption and (b) NH3 adsorption. The iso-value is 0.01 
e/A3. 

from the weak physisorption of gas molecules on the graphene surfacejisl. l27||. Furthermore, 



the Bader analysis |36l| of the charge distribution is used to understand the nature of the inter- 
action between the gas molecules and the AGNRs, and to evaluate the induced effects on the 
molecules. The trend of calculated charge transfer (Table 1) can be understood on the basis 
of relative electron-withdrawing or -donating capability of the adsorbed molecular groups. 
From Table 1 (note: the positive value means a charge transfer from the adsorbed molecule 
to the AGNR), we can see that CO, NO, NO2, O2, and CO2 have electron- withdrawing 
capability, while NH3 is electron-donating functional molecule. It is well known that NO2 
and O2 are relatively strong electron-withdrawing molecules while NHois relatively strong 



27 



29l | . Comparing 



electron-donating molecule in other carbon-based materials 0, [s, Q, 
with gas adsorption on the graphene (GNR) surface {2^, the values of the charge transfer 
are much larger. This is consistent with our conclusion that the interaction between the gas 
molecules and the GNR edges is much stronger than that of surface. 

The calculated band structures and DOS of 10-AGNRs with molecule adsorption are 
shown in Fig. 3. Comparing with the DOS of AGNR (Fig. lb), the total DOS of the system 
and LDOS of the molecules show that these molecules modulate the electronic property 
of AGNRs in different manners: i) CO and NO molecules adsorption introduces impurity 
states in the band gap and the Fermi levels of two systems cross these states, as shown 
in Figs. 3a and 3b. Therefore, gas adsorption will decrease the original band gap, and 
probably have some influence on the optical properties of AGNRs. For CO adsorption, 
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there are two half-occupied states in the band gap, but we do not expect these impurity 
states can enhance the conductance of the system because these states are very locahzed 
and deep in the original band gap. It would be very difficult for charge carriers to transit 
between the valence (or conduction) band and impurity states at finite temperature, ii) 
LDOS analysis (Fig. 3c) shows that NO2 adsorption will introduce fully-occupied states 
which are strongly hybridized with the original "bulk" states in the valence band and these 
states are nonlocalized. It suggests that the interaction between NO2 molecules and dangling 
bonds of the ribbon is very strong, consistent with the calculated adsorption energy. The 
Fermi level is pinned in the top of valence band, which is the same as the case of the 
ribbon without molecule adsorption, so the system is still semiconducting, iii) Figs. 3d 
and 3e show the cases of O2 and CO2 adsorptions, respectively. From LDOS analysis we 
can see that the states contributed by CO2 (or O2) molecules are localized around the 
top of valence band and hybridize with the original valence band. Partial charge density 
analysis (Fig. 4a) shows that the states near the Fermi level are quite localized and mainly 
contributed by CO2 molecule and the carbon atoms of ribbon around the CO2 molecule. 
This suggests that the conductance of this system can not be enhanced notably. When 
the molecular doping concentration is low enough, these impurity states will become more 
localized on the gas molecules. But due to these half-occupied impurity states being near the 
top of valence band, the electrons of the valence band can transit into these states and the 
system will exhibit ^-type semiconducting behavior at finite temperature, iv) NH3 molecule 
adsorption induces unoccupied local states in the conduction band, and more importantly, 
the Fermi level is shifted into original conduction bands, resulting in n-type semiconducting 
behavior (Fig. 3f). Furthermore, partial charge density analysis (see Fig. 4b) shows that 
the states near the Fermi level are mainly contributed by the carbon atoms of the ribbon 
rather than NH3 molecules. The above results indicate a transition from semiconducting to 
conducting behavior after NH3 molecule adsorption. The stabilization of the Fermi level in 
the conduction ban d by the impurity resonant levels was also reported in semiconductors 
doped by Cr and Tl|37|, l38 |. 

Among all gas molecules considered, obviously NII3 molecule adsorption can greatly 
enhance the conductance of AGNRs, where the system will exhibit metallic behavior after 
NH3 adsorption. CO2 and O2, on the other hand, may enhance the conductance of GNRs 
to some extent at finite temperature since the CO2 and O2 adsorption will turn AGNRs 
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FIG. 5: The /-14)ias curves for the GNR sensor before and after the adsorption of NH3 and C02- 
The inset shows the schematics of such a GNR sensor, consisting of one 10-AGNR (detection region) 
and two metalhc 7-ZGNRs leads. The gas molecules can be adsorbed around the DB defects of 
the GNR sensor. 

to ^-type semiconductor. Based on the analysis, we can expect that AGNRs may act as 
effective sensor to detect NH3 out of other gas molecules discussed above by measuring the 
change of conductance after the gas adsorption. What is more interesting, the adsorption 
energy of NH3 molecule is only -0.18 eV (Table I), so NH3 molecules can be desorbed at 
higher temperature. This implies that GNR senors could be recycled more than once. 

A design of a GNR-based junction (GNR sensor) to detect NH3 is given in the inset of 
Fig. 5 as an example. It contains a 8.60 nm long 10-AGNR with one DB defect per edge 
as the detection region and two semi-infinite metallic 7-ZGNRs as the leads. Gas molecules 
can be adsorbed on the DB defect sites, and the conductance is measured by applying a 
bias voltage through the junction. The DB defect concentration in this model is about 
0.011/A per edge, which is practical in experiments. We calculate a series of current versus 
bias voltage (/-Vbias) curves for such GNR junction with different gas molecule adsorption 
on the edges. For the sake of clarity, we only show / — Vbias curves for the AGNRs before 
and after NH3 and CO2 adsorption due to the fact that the currents induced by other gas 
molecules adsorption are almost zero (much smaller than the currents induced by CO2 or 
NH3 adsorption). As shown in Fig. 5, without gas molecule adsorption, the channel (GNR 
sensor) exhibits a semiconducting behavior, and the current is always zero even under a bias 
of 0.5 V. After NH3 adsorption, however, the current increases notably and I-Vuas curve is 
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nearly linear, corresponding to a metallic behavior of ohmic contact. CO2 adsorption can 
also increase the current, but its value is much smaller than that induced by NH3 adsorption. 
This phenomenon indicates that NH3 can be detected out of other gases by applying a bias 
voltage upon the GNR junction, which is consistent with our analysis based on electronic 
properties. 

IV. SUMMARY 

In summary, we have performed first-principles calculation to study the adsorption ge- 
ometries and electronic structure of graphene nanoribbons with gas molecule adsorption. We 
find that NH3 molecule adsorption can significantly influence the electronic and transport 
properties of AGNRs, while other gas molecules have little effect. Based on this charac- 
teristic, we demonstrate that an AGNR can be used to detect NH3 molecules out of many 
familiar gas molecules. Furthermore, our work also suggests an effective way to fabricate 
n-type (p-type) transistors by NH3 (CO2 or O2) adsorption on graphene nanoribbons. 
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